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Abstract 

Recently, a new method for calculating the Cherenkov field acting on a point¬ 
like electron bunch passing through longitudinally homogeneous structures lined 
with arbitrary slowdown layers was proposed, where the formalism was obtained 
though consideration of a general integral relation that allows calculation of the 
fields at the vicinity of a point-like bunch. It demonstrates that the Cherenkov 
field at the point of the short relativistic bunch does not depend on the waveguide 
system material and is a constant for any given transverse dimensions and cross- 
section shapes of waveguides. With this paper we present a strict derivation of the 
fields formulas valid at the cross-section of a bunch on the basis of a conformal 
mapping method. We generalize the results of the previous paper to the case of 
transversely distributed bunch by deriving a two-dimensional Green function at 
the cross-section of a bunch. Comparison of the results proving validity of the 
method is given in the appendix. 


1 Introduction 

Relativistic, high intensity and small emittance electron bunches are the basis of linear 
collider (ILC [^, CLIC [^) and EEL (LCLS, X-FEL, etc.) projects among many 
others. These bunches excite Cherenkov wakehelds when electrons pass through the ac¬ 
celerating structures or other longitudinally extended components of a beam line (pipes, 
collimators, bellows). Theoretical analysis of Cherekov radiation commonly considers a 
’’short bunch” approach It can be applied for various Cherenkov generation pa¬ 

rameters, where the moving bunch size is much less that the fundamental wavelength if 
the high frequency spectrum range is not under investigation. Recently a new theoretical 
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approach that can be used for obtaining direct analytical formulas for electromagnetic 
held components at the position of a point-like bunch was proposed [8 10 . It was 
demonstrated that the longitudinal part of the Lorentz force acting on the point-like 
bunch does not depend on the waveguide system material and is a constant for any 
given transverse dimensions and the cross-section shapes of waveguides. The equiva¬ 
lence and exact matching of the longitudinal electric helds on beams passing through 
various waveguide conhgurations was also subject to analysis. The proposed approach 
considers using an integral relation based on the cylindrical slow-wave structure model. 
For planar, square and other cross-section geometries, one can obtain a corresponding 
form factor multiplier by using a conformal mapping from the solution for the cylindrical 
case [8 -10 . However, the question about the transverse distributions of the helds and 
Lorentz force in the cross-section of a bunch was not examined in details. The method 
described in [sll^ gives the opportunity to calculate the longitudinal electric held 


and longitudinal derivative of the transverse part of the Lorentz force at the point of a 
bunch, which is also assumed to be point-like. In this paper we are going to extend the 
idea of [^|^ to a bunch which is a point only by longitudinal coordinate and has some 
arbitrary distribution by transverse coordinates and derive formulas that gives distribu¬ 
tion by the transverse coordinates of the longitudinal helds Ez, Hz and the transverse 
part of the Lorentz force in the cross-section of a bunch. This formulas are the Green 
function and therefore may be applied to the helds calculation of distributed bunches by 
transverse coordinates. We also show in Section]^ how the relativistic Gauss theorem [^ 
can be derived with the use of the proposed approach. 


2 Main formular derivation 


In this section, we are going to derive the main formula for the transverse distribution 
of the longitudinal helds Ez and Hz in the cross-section which includes a bunch. The 
method is based only on the assumption that a complete electromagnetic held in the 
cross-section of a bunch is equal to zero outside the vacuum channel. 

Let us consider the Maxwell system in GGS units inside the vacuum channel: 


V X E 

V X H 


V • E 

V H 


laH 

c dt 

dvr. 1 aE 

c dt ^ 

Anp, 

0 . 


( 1 ) 


We assume that the bunch is moving along the z-axis of the structure and the structure 
is longitudinally homogenous. In the ultra-relativistic limit the current that produces 
the bunch and it’s charge density could be written as: 


j = jz = cp, 

p = Q5{x-xo)5{y-yo)6{z-ct). (2) 

Let us introduce a new coordinate ( = ct — z and split the Maxwell system into two 
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parts: 


and 


[V X E]^ 
[V X H]^ 





(3) 


V_L X E_l 
X 


■ E_l 


V^-H^ 


dH, 

■ ac ’ 


dvrp + 


dvrp + 




’ 

dE, 

dC ’ 


dH, 

dC ■ 


(4) 


Here E -symbol indicates components orthogonal to the z-axis and V_l - is the nabla 
operator in the plane orthogonal to the z-axis. From the eqnations ([^ one can achieve: 

dH, _ dE, 

dx dy ' ^ ’ 

dH, _ _dE, 
dy dx 


Eqnation (|^ conld be considered as the Canchy-Riemann eqnation for some fnnction 
dehned in a complex plane. Let ns consider the waveguide vacuum channel cross-section 
which includes the bunch as a complex plane. Let us rewrite Q 


as: 


V_L ■ E^ -I- iV_L X Ej_ = dvrp -|- 


dE.. 


.OH. 


dC dC ’ 


( 6 ) 


Vi . Hx + 'iVi X Hi = * I 47r^)+ ^ 


We introduce complex functions: 


We introduce an operator [TT|: 


One can see that 


^ Ex 'f^Ey^ 
h — Hx E xHy. 


Vc 


A A 

dx dy 


(7) 

( 8 ) 


V*e = V_L ■ E_l + fV_L X E_l, 

V:h = Vj.-H^ + iV±xHj.. (9) 
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Here * - symbol means the complex conjugation. On the other hand one can rewrite 
as: 

d 


V. = 2- 


dx*' 

_ O ^ 

^ dx' 


( 10 ) 


here y = x + iy. From (|^ one can treat function L = as a plane vector 

held. According to 12 , as this held has no sources and vortexes, it could be found as a 
derivative of it’s complex potential: 

dgz dE, m. 


dx dC dC 

By combining ([^ , ([l0|) , d?]) and (11) one can derive: 

de dgz 

= 47rp + —. 

ox ox 


( 11 ) 


( 12 ) 


One can see that h = ie, thus we can consider only one equation for e. 

Let us consider a simply connected region D with the arbitrary smooth boundary 
and corresponding complex plane ca; with ([^ equation (12) could be written as: 


2^ = i7iQS(x - xo)S(y - yo)S(-Q + 


(13) 


According to the Riemann’s theorem 12 there exist a conformal mapping from the region 
D to a circle and a corresponding y - plane: y = x{^) with xi^o) = 0 (wq = xq + iyo). 
By introducing the substitution of variables in (13), we arrive at: 

dgz dx 


ay du 


dx 


du 


5{x')5{y')6{-C) + 


dx du 


(14) 


Here x' and y' are dehned as y = x' + iy'. Simplihcation of (p(4|) gives: 


dgz 

dx ' 


(15) 


From (15) we can see that the task of the derivation held e in the D region could 
be reduced to the task of the derivation e in a circle with the source of intensity 
AtiQ5{—C,) (^) placed in it’s center. Let us consider a circular contour F in the y 

cj=a;o 

- plane such that |y| = r < a (a -is the radius of a circle). According to 12 with (15), 
a contour integral over F could be written as: 


y) % 


(16) 


aj=a;o 


Now let us rely on the fact that in case of a cylindrical structure with the bunch driving 
along the center of a cylinder longitudinal components (and also their C, derivatives) of 
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the electric and magnetic field do not depend on the transverse coordinates. In 

this case, in x plane from 0 we have: 

9z = Cx, (17) 


where C is a constant. Dne to the symmetry of rotation in the y-plane the transverse 
field e conld be expressed as: 


e = e(|xl)^. 


(18) 


When placing (17) and (18) into (16) and taking into acconnt that |e| is a constant on 
r one can write: 


X 


c* 


\x\=r 


Ixl 


dx 


ndx-^ <p X*dx = 4:TriQ6{-C) ^ 


doo 


IxNi- 


(jJ=U}Q 


Evalnation of the integrals on the left side of (19) gives: 

2Q6{-C) dx 


e = 


du 


LJ=LJO 


C* 

H-r, 

2 


mnltiplication of e by and taking into acconnt that |x| = r gives: 


(19) 


( 20 ) 


2Qi{-C) 

X* du 


+ 


LJ=UJQ 


2 


( 21 ) 


As we are considering the cross-section that inclndes a bnnch and follow the idea in 
we conclnde that the field at the channel bonndary vanish dne to the phase speed of light 
in the snrronnding medinm being less than the speed of light in the channel, or we can 
say that the effective plane field gz cansed by the radiation prevents penetration of the 
free space field of the bnnch into the medinm. Taking this into acconnt, one can write 
down a condition for the nnknown constant C: 


e(|xl = a) = 0. 


Using this condition from (21) we find: 

4Q5(-C) f dx 


Thns: 


C = -- 


e = 2 QS (-0 ^ 

du 


du 


UJ=UJ0 


iJ=LJo 


X 


X a" 


From ( [II| , ( [I^ and ( [2^ we have: 


dE^ .dH^ _ AQ5{-C,)dx[dx 

sc “ sc “ 


a? du V du 


(jJ=!jjQ 


( 22 ) 


(23) 


(24) 


(25) 
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Finally at the point C = 0 (the cross-section that includes the bunch) we have: 


du V du 


( 26 ) 


UJ=iJQ 


if x{^) = f{^) ■ function that gives the mapping of D -region on a circle we can write 
result in a more compact form: 


(Jj 

= -^Im[/'(a;)*/'(cuo)]. 

(Jj 


(27) 


Formula (27) gives values of the longitudinal held components and in the cross¬ 


section of a bunch which is moving along the vacuum channel of arbitrary shape. The 
formula is valid at the point of and away from the bunch and is the two-dimensional 


Green function by the transverse coordinates. By integrating (27) over uq with the 


transverse distribution of a bunch p±{uJo), one can achieve helds for the bunch distributed 
in transverse coordinates. 


Another method for deriving formula (27) based on the Poynting’s theorem and more 


physical approach can be found in Appendix A.l 


3 Transverse part of the Lorentz force 

We consider a part of the Maxwell system ([^ and write it down in the expanded form: 


dE, dE„ 


OH, 


+ 

dC 

dE, 

dHy 

ac 

dx 

dC 

dH, 

dHy 

dEx 

dy 

^ dC- 

dC, 

dHx 

1 

CO 

1 

dEy 


d( dx d( 


By combining (28),(29) and (30),(31) we have: 


2—(F 

dC dy 


2^ (E, + H,) = 


dH, 

dx 


dE, 

dx 

dE,, 


dy 


(28) 

(29) 

(30) 

(31) 


(32) 


Using the dehnition of the Lorentz force acting on a bunch we introduce a complex 
function E\ : 


E± — q [Ex — Hy -\- i{Ey Hx )], 


(33) 
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where g - is the charge of the test particle. Let us rewrite (32) with (33) as: 


2 ^^ {Fx + iFy) — q 


dH, dE, 


dy dx 


— i 


dH, dE, 

+ 


dx dy 


In terms of the derivative (10) hj uj = x + iy 

dF* d 


the equation above simplihes to: 




In a vicinity of C = 0 with (25) we have 


F, = 


4gQC f 

dio'^ 


dx 

doj 


( 34 ) 


(35) 


(36) 


U}=U)Q 


If x{^) = f{^) ■ function that gives the mapping of the D -region on a circle we can 
write down the result in a more compact form: 


F, = 


4gQC 


nujTf'M. 


(37) 


Formula (37) gives values of the transverse part of the Lorentz force acting on a bunch 


which is moving along the vacuum channel of arbitrary cross-section shape. By transverse 
coordinates, the formula is valid at the point of and away from the bunch and is the two- 
dimensional Green function. By the longitudinal coordinate, the formula is valid in a 


vicinity of a point C = 0- By integrating (37) over uq with transverse distribution of a 


bunch p±{uJo) one can achieve the Lorentz force for the bunch distributed in transverse 
coordinates. 


4 Relativistic Gauss theorem 


In this section we are going to show how one can derive the relativistic Gauss theorem 
introduced in with the use of a formalism demonstrated in the Section 

Let us consider equation (13) and integrate it over the vacuum gap cross-section: 


df* f do* 

dxdy = AttQ5{—C,) + f 


duj 


doj* 


(38) 


Svac Svac 

Let us evaluate an integral on the left side according to the Green’s theorem: 

/■ de* / 

-dxdy = —i/2 ® e^dcn. 


duj 


(39) 


here F' - is the boundary of the vacuum channel region S^ac- 

Taking into account condition (22) (e=0 on F') with ([IT| we immediately have: 


^ I / E,dS + i I H,dS\ =-47 tQ6{-C). 

Svac 


( 40 ) 
























In a vicinity of the the point (^ = 0 we have: 

[ E,dS + i [ H,dS =-A7iQe{-C), 


( 41 ) 


here 9{—() is a Heaviside theta-function. 


Expression (41) is a slightly generalized form of the relativistic Gauss theorem formu¬ 
lated in However, equation (41) states that an integral over the vacuum channel 


from the magnetic held Hz is zero in the cross-section of the bunch and in a vicinity of 
the the point C = 0- This result was not mentioned in 1^. 


5 Conclusion 

In this paper we have derived general formula for the longitudinal held components 
Ez and Hz (27) in the cross-section of a bunch and a vicinity of the point C = 0. 


The formula expresses the held though the derivative of the function that conformally 
remaps the channel of the form one is interested in on a circle . It is worth to mention 
that the mapping function itself does not appear in the formula and one needs to know 
only it’s derivative. It simplihes the calculations even more as, for example, in case 
of the Christohel-Schwarz integral helds expressed in terms of the function under the 
Christohel-Schwarz integral. 


Now, if one take the point of the bunch (which corresponds to the a; = cjq in (27)), 
we arrive at the result of the paper [^: 

E,(0) = -^|/'Mt = -^|J|, (42) 

here \J\ is the determinants of the Jacobi matrix at the point of the bunch. 


We also derived a formula for the transverse Lorentz force (37), which could be used 
for the bunches with non-point transverse distribution. 
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Appendix A 

A.l Derivation based on the Poynting’s theorem 

Let us start from the Poynting’s theorem 


Here 


S=-[ExH], 
ti = T (ED + BH). 

OTT 


(Al) 

(A2) 
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Using the fact that the square of the transverse part of the Lorentz force divided by the 
charge = F±/Q inside the vacuum channel could be expressed as: 

h" = (A3) 


and introducing a new coordinate C, = ct — z one can rewrite (Al) for the vacuum channel 
as follows: 


d(El + Hi) dfl te 

2SC + W 

Taking into account the fact that the bunch is moving along the z -axis and 

j = Qc6{x - xo)6{y - yo)S{-C), 


equation (A4) takes the following form: 


■ Sx = -inQEJix - x„)Hy - !/„)i(-C). 


Now let us integrate (A6) over the cross-section which includes the bunch: 


d{El + Hl)_ ^ dfl 


dS ^ - f S±ndl = 


2dC 2dC 
=-4:7rQ6{-C) [ EJ{x -XQ)d{y -yo)dS. 


(A4) 


(AS) 


(A6) 


(A7) 


Taking into account that energy flux in the cross-section of the bunch is collinear to the 
beam direction, the orthogonal component of the Ponynting vector should be equal to 
zero as a consequence of the transverse Lorentz force in this cross-section being zero. By 
integrating over ( from — oo to zero, one can derive: 

J [E'^ + H'^'jdS =—AttQ j EzS{x — xo)5{y — yo)dS. (A8) 

Here integrals are taken over the vacuum gap cross-section. Let us also take into account 
the consequence of the Maxwell system: the equation for longitudinal components of the 
electric and magnetic field in an ultra-relativistic case inside the channel: 


^±E, = 0, (A9) 

A±Hz = 0 . 


Here Aj^ - is the transverse part of the Laplace operator. From equations (A9) one can 
conclude that E^ and are harmonic functions inside the vacuum channel and fulhl the 
Cauchy-Riemann equations (|^. This means that one can introduce a complex analytical 
function inside the vacuum channel 


w = Ez — iHz- 


(AlO) 
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The modulus square of this function according to (A8) will give energy density inside 
the cross-section: 


\w\^ = E^^+Hl (All) 

Let us consider the case of a cylindrical structure when the bunch is driving along 
the axis of the cylinder [^[^: 

E, = -\ (A12) 

Let us consider a conformal mapping of the arbitrary cross-section (w-plane )on the circle 
(x-plane): 


X = /(w)- 


(A13) 


In this case formula (A8) in the coordinates of the y - plane could be written as: 


[El + Hl]\f'{^)fdS^ = in / EJlrix^x^lWflinlldS- 


(A14) 


An integral on the right side is equal to 

4nQ j E,5[r{x^,X2)]\f\uj)\^dS^ = AT^QE,\f\uj,)\\ (A15) 

Here Uq - is the point in the u - plane that corresponds to the center of the circle in the 
y - plane. On the other hand in the y - plane we have: 

A'kQ j Ez6[r{xi,X2)\dS^ = A'kQE;, = j [El + E[l'\dS^, (A16) 


thus 


Et + Ht 

Z ' Z 


dS^ = |/'(a.„)|" / [El + Hi] |/'(a.)|"dS“ 


(A17) 


Taking into account that E^ held in the y - plane is dehned by (A12) and = 0 , one 
can write: 


Ef + Hi 


dS‘^ = 


4g2 


\fM\^ / \f{u:)fdS‘^ 


From (All) we conclude that 


w 


2Q 




From we know that 


M;((:no) = E^ 


2Q 


i/'K)r. 


(A18) 


(A19) 


(A20) 
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This means that there are only two possible variants of the function w. 



(A21) 


(A22) 


Equations (A21) and (A22) differ by the sign of a longitudinal component of the magnetic 
held. Taking into account ([^ and the fact that /(ca) is an analytic function and satishes 
the Cauchy-Riemann equations, we conclude: 


E, = -^Re[/'(a;)7'M], 

CL 


One can see that the result (A23) is equal to the result (27) of a Section]^ 


(A23) 

(A24) 


Appendix B 


In this appendix we present results of the transverse structure calculations of the longi- 



nique to derive a held in case metal walls are present in the cross-section of a bunch. 
Rectangular waveguide with side metal walls is considered as an example and formula 
for the longitudinal electric held Ez at the point of a bunch in case the bunch is driving 
in the center of the structure derived. 


B.l Transverse distribution of the field in a cylindrical waveg¬ 
uide 

Conformal mapping of a circle with the radius a on to a circle with radius a such that 
the point ujo = tq corresponds to the center of the second circle is given by 


Derivative is given by: 


/(^) 


= a 


CO -rp 
— ujTp 


n^) 


— Tq) 

(a^ — Tpuy' 


Assuming u = rexp[i0] one can achieve: 


(B25) 


(B26) 
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1 2 3 4 5 6 


Figure Bl: Angular dependence of in case of the displaced bunch and r = tq. Red 
dots are for the direct simulation using the mode decomposition method 13,14 and a 
solid line is for the formula (B29). 



Figure B2: Angular dependence of in case of the displaced bunch and r = a (on the 
channel boundary). Red dots are for the direct simulation using the mode decomposition 
method [T^ITd] and a solid line is for the formula (B29). 
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1 2 3 4 5 6 


Figure B3: Angular dependence of in case of the displaced bunch and r = a/10. Red 
dots are for the direct simulation using the mode decomposition method 13 ,^ and a 
solid line is for the formula (B29). 


Table Bl: Parameters for the cylindrical structure. 


a 

b 

£ 

ro 

Q 

0.3 cm 

0.32 cm 

5.7 

0.15 cm 

1 nC 


f'M = 


„2 ’ 


(B27) 


and 


Re|/'(i..)-| = 


From (A23) we have: 


E.= 


a?{a? — 'rQ)(a^ — 2a?rrQ cos(0) + r^rg cos(20) 
(a^ + r^rg — 2a^rrQ cos(0))^ 


2Q a® — 2a®rro cos{(j)) + r^rgo"^ cos(2(/) 
(a^ + r^rg — 2a?rrQ cos(0))^ 


(B28) 


(B29) 


Table jBT] contains a list of parameters for the cylindrical structure. Results are de¬ 
picted in Figure jBTj Figure jB^ and Figure jB^ Figures show angular dependence of the 
longitudinal electric field E^ for three different values of the test charge radial coordi¬ 


nates. We see full agreement with the mode decomposition method 14 
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Q _I_I_I_I_I_I_I_1_I_ 

0.5 0.9 1.3 1.7 2.1 2.5 2.9 3.3 3.7 4,1 4.5 

cm 


Figure B4: Y-dependence of Ez in case the test charge is located a.t x = 0.4 cm from the 
center. Red dots are for the direct simulation using the mode decomposition method 
and a solid line is for the formula (B34). 


B.2 Transverse distribution of the Ez field in a planar waveg¬ 
uide 


Conformal mapping of a strip with the half-high a on a circle with the radius a such that 
the point cjq = 0 corresponds to the center of the circle is given by: 


Derivative is given by: 


f{oj) = a tan 


'ITU' 

- 4a - 



'TTOJ' 

- 4a - 


2 


(B30) 


(B31) 


Table B2: Parameters for the planar structure. 


a 

b 

e 

ro 

Q 

0.5 cm 

0.502 cm 

5.7 

0 cm 

1 nC 


Assuming u = x + iy, one can achieve: 


and 


Re|/'(‘^)*l = 


fM = 1 

(B32) 

TT 1 + cos [f ] cosh [^] 

^ (cos [f] + cosh [^]Y 

(B33) 
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cm 


Figure B5: X-dependence of Ez in case the test charge is located sA y = 0.7 cm from the 
center. Red dots are for the direct simulation using the mode decomposition method 
and a solid line is for the formula (B34). 


From (A23) we have: 


Ez = 


2Qn‘^ 1 + cos [f] cosh [|^] 

^ (cos [f] + cosh [If])' 


(B34) 


Table jB^ contains a list of parameters for the planar structure. Results are depicted 
on Figure jB4| and Figure jB5j Figures show coordinate dependence of the longitudinal 
electric field Ez on the x and y coordinates. We see full agreement with the mode 
decomposition method (Isj. 


B.3 Longitudinal electric field Ez on a bunch in a rectangular 
structure with side metal walls 

We consider a rectangular structure which differs from the planar by two parallel perfect 
conducting walls. In case the field distribution in the vacuum gap is known for the 
planar case (B34), contribution of the metal walls could be included by introducing 
mirror charges. If the field for the planar case is Ef{x,y) and metal walls are placed 
sA y = -Ew/2, a full field could be found as a superposition of the real charge field and 
imaginary charges held. In case the charge is placed in the centre of the structure, a full 
held could be found as: 

OO OO 

E:‘‘(x,y) = Et(x,y) + Y,(-irEi(x,y - nw)+ Y^(-irE?^(x,y + nw). (B35) 


n=l 


n=l 


The logitudinal electric held Ez at the point of a bunch could be found as: 


^rec 


(0, 0) = EPj{0, 0) + 2 (0, nw). 


(B36) 


n=l 
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Figure B6: Normalized longitudinal electric field on a bunch in the rectangular waveg¬ 
uide; dependence from the waveguide relative width a/w. Red line - dependence of for¬ 
mula (B37) on a/w normalized by — black dots - same dependence achieved using 
the mode decomposition method 15 and also normalized by — 


By placing (B34) into (B36), one can achieve: 




(- 1 )' 


i 1 + cosh (Sg!!) 


(B37) 


Figure |B^ shows dependence of the normalized longitudinal electric field at the point 
of the bunch on a relative width of the waveguide - a/w. One can see full agreement 


with the mode decomposition method 15 
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